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Abstract 



We point out that perturbative evolution in QCD at three loops generates a strange- 
antistrange asymmetry s(x) — s(x) in the nucleon's sea just from the fact that the 
nucleon has non- vanishing up and down quark valence densities. The recently 
computed three-loop splitting functions allow for an estimate of this effect. We find 
that a fairly sizable asymmetry may be generated. Results for analogous asymmetries 
in the heavy-quark sector are also presented. 
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1 Introduction 



Strange quarks and antiquarks play a fundamental role in the structure of the nucleon [1J. Among 
the various strangeness-related properties of the nucleon, the strange "asymmetry", s(x) — s(x), 
in the number densities of strange quarks and antiquarks, x being the light-cone momentum 
fraction they carry, is of particular interest. Since the nucleon does not carry any strangeness 
quantum number, the integral of the asymmetry over all values of x has to vanish: 

(s-s) = f 1 dx \s(x) - s(x)] = . (1) 
Jo 

However, there is no symmetry that would prevent the x dependences of functions s(x) and s(x) 
from being different. One therefore can expect s(x) ^ s(x), in general. 

To understand and quantify the strangeness asymmetry in the nucleon is interesting in various 
different contexts. Some models of nucleon structure make predictions EJ |4j 13 [U for 
s(x) — s(x), and their confrontation with experimental measurements may perhaps give us further 
insight into non-perturbative dynamics of the strong interactions. For example, within the meson 
cloud model 013113, contributions to the strange sea arise from fluctuations in the proton wave 
function to AK and Y>K states. The A (which contains the s quark) will primarily carry the larger 
fraction of the nucleon momentum, and one thus expects s(x) to be larger than s(x) at large x, 
implying the opposite behavior at small x. However, a number of subtleties in this picture have 
been pointed out and discussed in the literature [:9|. Light-cone models [5|, on the other hand, 
generically predict s(x) — s(x) < at large x. The various models have in common that they 
predict a fairly small value of the second moment of the strange-antistrange distribution, 

(x(s — s)) = / dxx \s(x) — s(x)] , (2) 
Jo 

usually \ (x(s — s))\ ~ 10 -4 . 

As was emphasized in ifTDlfTTIl . the question concerning the strange asymmetry in the nucleon 
becomes particularly relevant in view of the "anomaly" seen by the NuTeV collaboration in their 
measurement of the Weinberg angle in deeply-inelastic neutrino scattering. The NuTeV result 

G2, 

sin 2 #w|NuTeV = 0.2277 ± 0.0013 (stat .) ± 0.0009 (sys .), (3) 

deviates around three standard deviations from the commonly accepted value sin 2 # w = 0.2228 ± 
0.0004 lfT3l . This large difference could be at least partly explained HOjJTTJ by a positive value 
of {x(s — s)). Typically, a value (x(s — s)) 0.005 would be required if one wanted to attribute 
the NuTeV anomaly to the strange asymmetry alone. Note that, if (x(s — s)) is indeed positive 
and if one assumes that s(x) — s(x) has only one node, then the vanishing of (s — s) implies that 
s(x) — s(x) is positive at large x. 

The NuTeV Collaboration has itself determined the second moment (x(s — s)) from a lowest- 
order QCD analysis of neutrino dimuon data [ 14 1 and finds a negative value 0*511 : 

(x(s-s)) = -0.0027 ±0.0013 . (4) 

Such a value increases the discrepancy in sin 2 6* w to a 3.7a effect. The second moment had also 
been investigated in "global analyses" of unpolarized parton distributions. Ref. lfT6l reported an 
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improvement in the global analysis if the asymmetry s(x) — s(x) is positive at high x. They 
found (x(s -s)) = 0.002 ± 0.0028 at Q 2 = 20 GeV 2 from their best fit. A recent update of this 
analysis ifTTl . however, reduces the asymmetry significantly. The most recent global QCD fit lfl~8l 
finds a large uncertainty for the strange asymmetry and quotes a range —0.001 < (x(s — s)) < 
0.004. 

The discussion reported so far regards strange-anti strange asymmetries that do not depend 
on the hard- scattering scale Q at which the nucleon is probed. In this letter, we point out that 
perturbative QCD definitely predicts a non-vanishing and Q-dependent value of the strange- 
antistrange asymmetry. We will show that non-singlet DGLAP evolution of the parton densities 
at three loops generates a strange asymmetry even if it is not present at the input scale for the 
evolution. Thus, we can provide a prediction for the strange asymmetry s(x) — s(x) based solely 
on perturbative QCD. The effect arises because at that order of perturbation theory the probability 
of a splitting q — > q' becomes different from that of q — > q', and because the nucleon has u and d 
valence densities. The three-loop splitting functions have very recently been published lfT9ll20ll . 
among them the splitting function needed for our perturbative estimate of s(x) — s(x). To begin 
with, we write down the evolution equations and determine the solution for the generated strange 
asymmetry. We then present some numerical results for the strange asymmetry and extend the 
analysis to the heavy-quark sector. 



2 Non-singlet evolution equations and their solutions 

The parton distributions f a (x, Q 2 ) (a = q^ q\, g) of the nucleon evolve according to the evolution 
equation 

= E jf 7 ■<*«">) W 

where P ab is the function describing the splitting b — > a. The splitting functions are perturbative; 
their perturbative series starts at 0(as): 



S)^ ,+ (S)^ ,+ (S)^ ,+ °^)' < 6) 

Keeping just the first term yields the leading order (LO) evolution. Improving the approximation 
by taking into account also the second, or the second and third, terms corresponds to next-to- 
leading order (NLO) and next-to-next-to-leading order (NNLO) evolution, respectively. 

It is convenient to introduce Mellin moments of the parton densities and splitting functions, 

/f(Q 2 )= P dxx^fa&Q 2 ), (7) 







so that momentum-fraction convolutions as in Eq. © reduce to true products and the evolution 
equation becomes simply 

d -§^ = E P& («s(Q 2 )) f b N (Q 2 ) ■ (8) 

In the following we drop the dependence on the moment index iV for simplicity. 
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Since s(x) — s(x) is a flavor non-singlet (NS) quantity, we only need to consider the NS 
sector of evolution. In the following we write the evolution kernels P ab by adopting the notation 
of Ref. |[T9ll . Owing to charge conjugation invariance and and flavour symmetry of QCD, one has 
(see e.g. Ref. ED) 

p p x tdV _i_ p5 

r giQk ~ r QiQk ~ U ik* qq T f qq 

PqiQk = Pq~i<lk = ^ikPqq + Pqq ■ (9) 

The splitting functions P^ q and P^ thus describe splittings in which the flavor of the quark 
changes. As will become clear below, the effect we wish to investigate originates from the fact 
that P^ ^ P ? 5 f starting from NNLO ED El- 

In the flavour NS sector the evolution equations © are diagonalized by properly introducing 
NS combinations of parton densities. Up to NLO it is sufficient to consider two NS combina- 
tions. Owing to the difference between P^ q and PjL starting from NNLO it is necessary ll22ll to 
introduce the following three NS combinations of parton densities: 

N f N f 
i=l 7V / j=l 

where Nf is the number of flavors. Each of these evolves as 

^™ = i«Wfl")) , (A=V,±), (11) 
where the evolution kernels are 

P (V) = K - ?l + Nf {Pq S q - P£) , P (±) = Pi ± Pi ■ d2) 

The equations have the solutions 

f (A) (Q 2 ) = U (A) (Q,Q ) f (A) (Ql) , (13) 

where f^ A \Ql) is the parton density at the starting scale Q and the evolution operator is 
given by 

U^(Q,Qo)=ex P \[ Q 2 2 ^P ( - A \a s (q"))) . (14) 
Using Eq. (THl) with A = — and A = V, we have 



(15) 

We remind the reader that the parton distributions as well as the evolution operators depend on 
the Mellin moment N. 



3 Strange-quark asymmetry 

Equation < fl3T > is the basic result in our discussion of flavour asymmetries. The key point is that 
Eq. (fT5l) implies that, in the region of Q 2 where QCD perturbation theory is applicable, the flavour 



3 






(a) " (b) 

Figure 1 : Example of (a) virtual and (b) real diagrams contributing to -P^ 5 



asymmetries (f q . — f qi ) (x, Q 2 ) must necessarily be different from zero (they can vanish, at most, 
at a single value of Q). This is a definite, though qualitative, prediction of perturbative QCD. 

In the following we simplify the notation, using f qi = qi and f q . = q h and we consider in 
detail the strange-quark asymmetry, s — s. Equation (fT5l gives 



S) (Q 2 ) = U(-\Q,Q ) 



(16) 



At LO and NLO, = U^~\ and thus any strange-quark asymmetry can only be produced 
by a corresponding asymmetry at the input scale Q of the evolution. Starting from NNLO, the 
degeneracy of P^ and P^ is removed: 

P (V) ~ P ( - ] = N f (P q S q - P$ = (g) 3 PS 8 + O(at) , (17) 

where P$ s has recently been presented in Ref. |fl9"l . It comes with the color structure d abc d a b c , 
which is also new at this order. In a physical gauge, the Feynman diagrams contributing to 
P( 2 J S are of the "light-by-light" scattering type, three gluons connecting the two different quark 
lines. Figure [T] shows some examples of (interferences of) (a) virtual and (b) real diagrams that 
generate the asymmetry in the evolution of quarks and antiquarks. The virtual part (e.g. Fig.[TJa)) 
has separately been studied [ 23 1 in the context of its contribution to the one-loop triple collinear 
splitting function. When the evolved quark qj is replaced by the antiquark <jj, the abelian-like 
part of the diagrams in Fig. [T] changes sign, because of the charge asymmetry produced by the 
exchange of three gluons (vector bosons) in the t-channel. This effect occurs both in QCD and 
QED, and it is a genuine quantum (due to interferences and loop contributions) phenomenon. 

The full expressions for Pj 2 ^ in N and x space may be found in Eqs. (3.9) and (4.11), 
respectively, of Ref. j\9\ ; also a simple approximation of the function in x space is provided 
which is sufficiently accurate for our purposes: 

P { n) S ( x ) ~ N f ([L^-WS^x^ 1 - 7.208 x) + 151.49 + 44.51 x- 43.12 x 2 + 4.82 x 3 ] [1 - x] 

+ L Li [-173.1 + 46.18 L ] + 178.04 L + 6.892 L 2 + 40/27 [L\- 2 Ljj]) , (18) 

where L = lnx, L\ = ln(l — x). The Mellin moments of this expression are straightforwardly 
derived and involve only simple harmonic sums [ 19 1. The first moment of the above parameteri- 
zation is zero to high accuracy. 
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On account of Eqs. dToll-dTSl). even if (s — s)(Ql) = 0, a non-vanishing strange-quark asym- 
metry is produced just by the perturbative QCD evolution. Here it is crucial that the total valence 
density of the nucleon, f^ v \ is non- vanishing due to the up and down valence quarks. Using 
Eqs. (fl4h and (flTl) we have, in moment space, 



u (v KQ,Qq) 

U(-\Q,Q ) 




P® s + (9(N 3 LO) 



47T 



where 



1 



12vr 



(l\C A - 2Nf 



O(N 3 L0) , (19) 



(20) 



and we have used the renormalization group equation 

dliaas(q 2 ) 



d\nQ 2 



P(a s (q 2 )) , 
-b a s - h al - b 2 a| 



0(ai 



(21) 
(22) 



Note that despite being a NNLO effect, the perturbative generation of s — s is a leading ef- 
fect since it first occurs at this order. This explains why only the coefficient b needs to be 
taken into account in Eq. (fl"9l) . Under the assumption (s — s)(Qq) = 0, and neglecting heavy 
quark (valence)contributions and threshold effects, the solution for the evolution equation for the 
strange-quark asymmetry reads to NNLO: 



>N 



(Q 2 ) 



P 



(2)5 



ns.N 



lTrb N f 



D (2)S 
ns.N 



8nbnN 




4tt J 
4tt / 



U i N \Q,Qo)f N v \Ql) 

(u^+d^) N m , (23) 



where we have restored the Mellin moment index N. In the second line we have used that 
(Q, Qo) simply evolves the valence input, = u — u and dP^ = d — d, to the scale Q at 
LO accuracy. Note that, assuming isospin symmetry, the sum of valence distributions is the same 
in the proton and the neutron and, consequently, also the perturbative strange asymmetry. 



4 Numerical estimates 

From Eq. (l23l) we can straightforwardly obtain predictions for [s—s] (x, Q 2 ) by a numerical Mellin 
inversion, once we have chosen an initial scale and input valence densities. For our estimates we 
employ the low input scale Q = 0.51 GeV and the u, d valence densities of the LO "radiative" 
parton model analysis of Ref. ll24ll . Threshold effects at Q = m c = 1.4 GeV and Q = m b = 4.5 
GeV are taken into account by the full implementation of Eq. (fl~6b . Since we are considering 
a leading effect, the LO approximations are appropriate. The value for the initial scale is of 
course crucial for our results; the lower the scale, the larger will the perturbatively generated 
strange asymmetry be at a given higher scale Q. Our choice of a rather low input scale may 
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be regarded as yielding the largest possible perturbative strange asymmetry. Whether or not it 
is indeed correct to assume that the nucleon is symmetric in s and s at a low scale is an open 
question; however, our motivation is to explore the new effect provided by three-loop evolution. 
We note that our input assumption [s — s] (x, Ql) is consistent with the input in Ref. ll24l . where 
actually s(x, Q%) = s(x,Ql) = was assumed, resulting in a purely generated (symmetric) 
strange distribution which agrees reasonably well with the ones obtained in other global analyses 
of PDFs 1 . 

Fig. |2 a) shows [s—s] (x, Q 2 ) as a function of x, for three different scales, Q 2 = 2, 10, 100 GeV 2 . 
For comparison, Fig. |2b) shows the ratio of [s — s](x, Q 2 ) to the MRST[25] strange density 
s(x, Q 2 ). One can see that the generated asymmetry is not negligible and turns out to be positive 
at small x and negative at large x. Since the distribution has a vanishing first moment and only 
one node, a negative second moment results: 

(x(s - s)) ss -5 x 10~ 4 (Q 2 = 20 GeV 2 ) . (24) 

This value depends fairly little on Q 2 once Q 2 > 1 GeV 2 ; it then very gently decreases at large 
Q 2 . As expected for a NNLO effect, it is quite small, somewhat smaller than the NuTeV value 
quoted in Eq. ©. We also note that our value lies in the band for the second moment derived 
from a global analysis in Ref. lfT8ll . 

Let us now try to put the size of the perturbatively generated strange asymmetry into a better 
perspective. As we discussed above, the effect becomes possible for the first time at NNLO, 
where -P^ 5 7^ 0. This is reminiscent of a well-known effect that first arises in NLO evolution, 
namely the perturbative generation of [u — d] (x) ^ from evolution, due to P^P ^ at NLO 
[26|. Interestingly, despite being a NNLO effect, we found our [s — s}(x, Q 2 ) to be larger than 
the NLO perturbative [u — d\(x, Q 2 ) in most of the x region, in particular at small x where the 
splitting function P^ s is singular as In 4 x (see Eq. (fT8t ) [fl~8l . Also, for the perturbative u — d 
the difference [u^ — d^](x, Qq) of input valence densities determines the boundary condition, 
whereas for s — s it is their sum, according to Eq. d23t . From this point of view, the perturbatively 
generated [s — s] (x, Q 2 ) can actually be considered as quite large. Of course, as it is well known, a 
much larger u—d asymmetry than the perturbatively predicted one has been measured 1271 12811291 
|30l|3T), which implies that non-perturbative effects outweigh the asymmetry from perturbative 
evolution. It is clearly possible that also in the case of s — s non-perturbative effects dominate. 
It is also worth pointing out that the uncertainties in the perturbative strange asymmetry itself are 
difficult to quantify since it is effectively a LO effect. On the other hand, as we mentioned in the 
introduction, models of nucleon structure generally predict a very small strange asymmetry, the 
second moment usually being several times smaller than ours in Eq. (l24l . Therefore, at the very 
least, we expect perturbative evolution to play a significant role in relating model predictions at 
the (low) model scale to s — s at scales relevant for comparison to experimental data. 

We also note that the large-x behavior of our perturbatively generated strange asymmetry is 
driven by that of the evolved valence densities and of the splitting function P^ s (x). As x — ► 1, 
the latter behaves as 

Pif^x) ~ (1 - x) ln(l -x) + 0(1 - x) , (25) 
which because of the convolution with the evolved valence densities implies that 

[s - s](x, Q 2 ) ~ (1 - x) 2 ln(l - x) \u {v) + d (y) l (x, Q 2 ) . (26) 



^This agreement supports the assumption that the asymmetry should also vanish at a low Qq 
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Figure 2: (a) Strange asymmetry in the nucleon generated from NNLO QCD evolution for Q 2 = 
2, 10, 100 GeV 2 and (b) the corresponding ratio to the LO strange distribution from Ref. K51 



This may well be the dominant behavior at high x, even in the presence of a non-perturbative 
input distribution for s — s. 

We finally note that our analysis may also be extended to predict the asymmetries for heavy 
flavors c and b. Here, the perturbative prediction may be more reliable since one will typically 
start the evolution from the mass of the heavy quark, which is in the perturbative region. As- 
suming that the charm and bottom densities vanish at the respective masses, we find the results 
shown in Fig.OJ The upper plot (a) compares the purely perturbatively generated charm and bot- 
tom densities to the results of the latest MRST LO global analysis ll25ll . The agreement found at 
scales far away from the threshold for heavy quark production, and in the relevant small x range, 
is a signal of the validity of the approach. One can expect that a similar situation holds for the 
asymmetry. The lower plot (b) corresponds to the ratio between the generated asymmetry and the 
corresponding heavy quark density. Note that for the last result we only assume that the heavy- 
flavor asymmetries vanish at the respective masses. This is a weaker assumption than that there 
be no initial heavy-quark distributions at all at these scales. As can be observed, the asymmetries 
are smaller than in the strange sector, mostly due to the larger initial scale at which the evolution 
begins. 



7 



(a) QP ertur ^ at ^ ve 




Charm / - 


7 Qmrst 




Bottom I: 






"lo- 2 


io ; 






Q 2 =100 GeV 2 - 


: (b) Q-Q 






: Qmrst 






3 


10 2 


io - 1 



X 

Figure 3: (a) Ratios between the purely pertubatively generated charm and bottom densities at 
Q 2 = 100 GeV 2 to the corresponding distribution from Ref. and (b) (Normalized) charm 
and bottom asymmetries in the nucleon generated from NNLO QCD evolution. 
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